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ABSTRACT 

We consider the problem of the thermal X-ray radiation from the hot polar cap of 
radio pulsars showing evidence of E x B subpulse drift in radio band. In our recent 
Paper I, using the partially screened gap (PSG) model of inner acceleration region 
we derived a simple relationship between the drift rate of subpulses observed in a 

■ radio-band and the thermal X-ray luminosity from polar caps heated by the back-flow 
I particle bombardment. This relationship can be tested for pulsars in which the so- 

fSJ , called carousel rotation time P4, reflecting the E x B plasma drift, and the thermal 

fN| ■ X-ray luminosity L x from the hot polar cap are known. To test the model we used 

only two available pulsars: PSRs B0943+10 and Bl 133+16. They both satisfied the 
model prediction, although due to low photon statistics the thermal component could 
not be firmly identified from the X-ray data. Nevertheless, these pulsars were at least 
cj . consistent with PSG pulsar model. 

O ^ In the present paper we consider two more pulsars: PSRs B0656+14 and B0628-28, 

whose data have recently become available. In PSR B0656+14 the thermal radiation 
from the hot polar cap was clearly detected, and PSR B0628-28 also seems to have 
such a component. 

I In all cases for which both P4 and L x are presently known, the PSG pulsar model 

seems to be fully confirmed. Other available models of inner acceleration region fail 
to explain the observed relationship between radio and X-ray data. The pure vacuum 
gap model predicts too high L x and too low P4 , while the space charge limited model 

■ predicts too low L x and the origin of the subpulse drift has no natural explanation. 
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1 INTRODUCTION 

Although almost 40 years have passed since the discovery of 
pulsars, the mechanism of their coherent radio emission is 
still not known. The theory of pulsating X-ray emission also 
demands further development. The puzzling phenomenon 
of drifting subpulses is widely regarded as a powerful tool 
for the investigation of the pulsar radiation mechanism. Re- 
cently, this phenomenon received a lot of attention, mostly 
owing to the newly developed techniques for the analysis of 
the pulsar radio emission fluctuations (Edwards & Stappers 
2002, 2003). Using these techniques, Weltevrede, Edwards, 
& Stappers (2006a, WES06 henceforth) presented the re- 
sults of the systematic, unbiased search for the drifting sub- 
pulses and/or phase stationary intensity modulations in sin- 
gle pulses of a large sample of pulsars. They found that the 
fraction of pulsars showing evidence of drifting subpulses 
is at least 55 % and concluded that the conditions for the 



drifting mechanism to work cannot be very different from 
the emission mechanism of radio pulsars. 

It is therefore likely that the drifting subpulse phe- 
nomenon originates from the so-called inner acceleration re- 
gion right above the polar cap, which powers the pulsar ra- 
diation. In the classical model of Ruderman & Sutherland 
(1975; RS75 henceforth) the subpulse-associated spark fila- 
ments of plasma circulate in the pure a Vacuum Gap (VG) 
around the magnetic axis due to the E x B plasma drift. 
This model is widely regarded as a natural and plausible 
explanation of the drifting subpulse phenomenon, at least 
qualitatively. On the quantitative level, this model predicts 
too high a drifting rate, or too short a period P4 (P3 in the 
nomenclature introduced by RS75), of the sparks' circula- 
tion around the polar cap, as compared with the observa- 
tions (e.g. Deshpande & Rankin, 1999; DR99 henceforth). 
Also, the predicted heating rate of the polar cap surface 
due to the spark-associated back-flow bombardment is too 
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high. The alternative model, namely the space charge lim- 
ited model (SCLF; e.g. Arons & Sharleman 1979), predicts 
too low a heating rate and has no natural explanation for 
the phenomenon of drifting subpulses (Zhang & Harding 
2000; Harding & Muslimov 2002). However, this model has 
an advantage over the VG model, namely it is free of the 
so-called binding energy problem, to avoid which the VG 
model requires an ad hoc assumption of the strong, non- 
dipolar surface magnetic field (for review and more detailed 
discussion see Gil & Melikidze 2002). 

Motivated by these observational discrepancies of the 
otherwise attractive VG model, Gil, Melikidze & Geppert 
(2003; GMG03 henceforth) developed further the idea of the 
inner acceleration region above the polar cap by including 
the partial screening caused by the thermionic ions flow from 
the surface heated by sparks. We call this kind of the inner 
acceleration region the "partially screened gap" (PSG hence- 
forth) Q Since the PSG potential drop is much lower than 
that in the RS75 model, the intrinsic drift rate Pa is com- 
patible with the observations. This is a consequence of the 
reduced potential drop, partially screened by the thermionic 
ion flow from the polar cap surface. In the pure vacuum 
RS75 gap, the heating of the polar cap is definitely too in- 
tense (e.g. Zhang, Harding & Muslimov 2000; Zhang, San- 
wal & Pavlov 2005; ZSP05 henceforth). On the other hand, 
the SCLF model predicts too low a heating rate as com- 
pared with observations (Zhang & Harding 2000; Harding 
& Muslimov 2002). Thus, by measuring the thermal X-ray 
luminosity from heated polar caps one can potentially re- 
veal the nature of the inner acceleration region in pulsars. 
This can also help to understand a mechanism of drifting 
subpulses, which appears to be a common phenomenon in 
radio pulsars. 

ZSP05 were the first who attempted to test different 
available models of the inner acceleration region in pulsars, 
using a concept of the polar cap heated by the back-flow 
particle bombardment. They observed the best studied drift- 
ing subpulse radio pulsar PSR B0943+10 with the XMM- 
Newton observatory and argued that the detected X-ray 
photons were consistent with PSG formed in the strong, 
non-dipolar magnetic field just above the surface of a very 
small and hot polar cap. Recently Gil, Melikidze & Zhang 
(2006 a,b; hereafter Paper I and II, respectively) developed 
a detailed model for the thermal X-ray emission from radio 
drifting pulsars. They applied their model to PSR B0943+10 
as well as to PSR B1133+16, which was observed in X-rays 
with Chandra observatory by Kargaltsev, Pavlov & Garmire 
(2006, KPG06 henceforth). These authors found that this 
case is also consistent with the thermal radiation from a 
small hot spot, much smaller than the canonical polar cap. 
PSR B1133+16 is almost a twin of PSR B0943+10 in terms 
of P and P values and, interestingly, both pulsars have very 
similar X-ray signatures, in agreement with the PSG model 
(see Table 1 and Fig. 1). 



1 Cheng & Ruderman (1980) were the first to consider the PSG 
model. However, they argued that even with partial screening in- 
cluded, the conditions above the polar cap are close to pure VG as 
in RS75. GMG03 demonstrated that the actual thermostatic self- 
regulation establishes the accelerating potential drop that may be 
as low as few percent of that of RS75 value. 



The PSG model can be tested if two observational quan- 
tities are known: (i) the circulational period P4 for drifting 
subpulses observed in the radio band (also called the pulsar 
carousel time), and (ii) the X-ray luminosity L x of thermal 
black-body (BB) radiation from the hot polar cap (see Eqs. 
2 and 3 below). The above mentioned observations of PSRs 
B0943+10 and B1133+16 are not decisive. Indeed, due to 
poor photon statistics, their spectra can be described by 
either a thermal model, a non-thermal model, or a combina- 
tion of the both. In any case, one can pose the upper limits 
for the thermal radiation from the hot polar cap from these 
data, so that the PSG model could be tested at least in the 
order of magnitude approximation. 

In this paper we include two more pulsars for which val- 
ues of both P4 and L x are currently known: PSRs B0656+14 
and B0628-28. The former case was a real breakthrough for 
our considerations and testing. Indeed, while in the other 
cases the character of the spectrum was not certain, in this 
pulsar (one of the Three Musketeers) the thermal radiation 
from the hot polar cap was clearly detected (De Luca et 
al. 2005). PSR B0628-28 was observed with Chandra and 
XMM-Newton observatories by Tepedelenlioglu & Ogelman 
(2005; hereafter TO05). We show that both pulsars comply 
the PSG model, increasing the number of pulsars that pass 
the model test expressed by Eqs. (2) and (3) from two to four. 
At the moment, PSRs B0943+10, B1133+16, B0656+14 and 
B0628-28 are the only pulsars for which both P4 and L x are 
known. It is important to show that all of them follow the 
theoretical prediction curve in Fig. 1. 



2 PSG MODEL OF THE INNER 
ACCELERATION REGION 

The charge depleted inner acceleration region above the po- 
lar cap results from the deviation of a local charge density 
p from the co-rotational charge density (Goldreich & Julian 
1969) pgj = ^ ■ B s /27rc « B s /cP. For isolated neutron 
stars one might expect the surface to consist mainly of iron 
formed at the neutron star's birth (e.g. Lai 2001). Therefore, 
the charge depletion above the polar cap can result from 
binding of the positive 26 Fe i° ns ( a t least partially) in the 
neutron star surface. If this is really possible (see Mendin & 
Lai 2006, and Paper II for details), then the positive charges 
cannot be supplied at the rate that would compensate the 
inertial outflow through the light cylinder. As a result, a sig- 
nificant part of the unipolar potential drop develops above 
the polar cap, which can accelerate charged particles to rela- 
tivistic energies and power the pulsar radiation mechanism. 

The ignition of cascading production of the electron- 
positron plasma is crucial for limitation of the growing po- 
tential drop across the gap. The accelerated positrons will 
leave the acceleration region, while the electrons bombard 
the polar cap surface, causing a thermal ejection of ions. 
This thermal ejection will cause partial screening of the ac- 
celeration potential drop AV corresponding to a shielding 
factor T) = 1 — pi I pgj (see GMG03 for details), where pi is 
the charge density of the ejected ions, AV = n(2iv / 'cP)B B h 2 
is the potential drop and h is the height of the acceleration 
region. The gap potential drop is completely screened when 
the total charge density p = pi + p+ reaches the co-rotational 
value pgj- 
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GMG03 argued that the actual potential drop AV 
should be thermostatically regulated and there should be es- 
tablished a quasi-equilibrium state, in which heating due to 
electron bombardment is balanced by cooling due to ther- 
mal radiation. The quasi-equilibrium condition is Q coo i = 
Qheat, where Q coo i = is the cooling power surface 

density by thermal radiation from the polar cap surface 
and Qheat = ym e c 3 n is the heating power surface den- 
sity due to back-flow bombardment, 7 = eAV/m e c 2 is the 
Lorentz factor, n = ncj — rii — r/ncj is the number den- 
sity of the back-flowing particles that deposit their kinetic 
energy at the polar cap surface, 77 is the shielding factor, 
rii is the charge number density of the thermionic ions and 
noj = PG.j/e = 1.4 x 10 11 &P°i 5 P"°' 5 cm" 3 is the corota- 
tional charge number density. It is straightforward to ob- 
tain an expression for the quasi-equilibrium surface temper- 
ature in the form T s = (6.2x 10 4 K)(P_i 5 /P) 1/ V /2 & 1/ V /2 , 
where the parameter b — B a jB a = A pc /Af, i describes the 
domination of the local actual surface magnetic field over 
the canonical dipolar component at the polar cap, and P_is 
is the normalized period derivative. Here A pc = 7rr 2 c and 
Abol = A p = TiTp is the actual (bolometric) emitting surface 
area, with r pc and r p being the canonical (RS75) and the 
actual polar cap radius, respectively. Since the typical polar 
cap temperature is T s ~ 10 6 K (Paper II), the actual value 
of b must be much larger than unity, as expected for the 
highly non-dipolar surface magnetic fields. 

The accelerating potential drop AV and the perpendic- 
ular (with respect of the magnetic field lines) electric field 
AE, which causes E x B drift, must be related to each other, 
and this relationship should be reflected in combined radio 
and X-ray data of pulsars showing drifting subpulses. This 
is basically a conal phenomenon (Rankin 1986), so we can 
restrict ourselves to the periphery of the polar cap, where 
these two potential drops are numerically equal to each 
other. Moreover, following the original "pillbox" method of 
RS75 we can argue that the tangent electric field is strong 
only at the polar cap boundary where AE = 0.5 AT/ /h = 
r](Tv/cP)B a h (see Appendix A in GMG03 for details). Due 
to the E x B drift the discharge plasma performs a slow cir- 
cumferential motion with velocity Vd = cAE/B s — r/nh/P. 
The time interval to make one full revolution around the 
polar cap boundary is Pt ~ 2nr p /vd- One then has 



If the plasma above the polar cap is fragmented into fil- 
aments (sparks), which determine the intensity structure 
of the instantaneous pulsar radio beam, then in principle, 
the circulational periodicity P4 can be measured/estimated 
from the pattern of the observed drifting subpulses (Desh- 
pande & Rankin 1999, Gil & Sendyk 2003). According to 
RS75, P4, — NPs, where N is the number of sparks con- 
tributing to the drifting subpulse pattern observed in a given 
pulsar and P3 is the primary drift periodicity (distance be- 
tween the observed subpulse drift bands). On the other hand 
N « 2nr p /2h = 7ra, where the complexity parameter can be 
estimated from the approximate formula a = 5P ( l'i^ P~ ' 64 
(Gil & Sendyk, 2000; GS00 henceforth). One has to real- 
ize that this approximation was derived under a specific as- 
sumption concerning the actual surface magnetic field (see 
discussion below equation (11) in GS00), and it can give 



misleading values of a for some untypical pulsars (see discus- 
sion in section 4) . However, using this concept we can write 
the shielding factor in the form 77 w (l/2ir)(P/Ps), which 
depends only on a relatively easy-to-measure primary drift 
periodicity P3. Note also that P4/P = a/{2rf). We show the 
values of the model parameters obtained from these equa- 
tions in Table 1. 

The X-ray thermal luminosity from the polar cap 
with a temperature T s is L x — aT*Tvr p — 1.2 x 
10 32 (P_i 5 /P 3 )(77/i/V P ) 2 erg/s, which can be compared 
with the spin-down power E = IQfl — 3.95/45 x 
10 31 P_i 5 /P 3 erg/s, where I = 7 4 5l0 45 g cm 2 is the neutron 
star moment of inertia0and I45 = llo' 2 2- Using equation JTJ 
we can derive the formula for thermal X-ray luminosity as 

L x = 2.5 x 10 31 (P_i 5 /P 3 )(P 4 /P)" 2 , (2) 



or in the simpler form representing the efficiency with re- 
spect to the spin-down power 




This equation is very useful for a direct comparison with 
the observations, since it contains only the observed quan- 
tities (although it is subject to a small uncertainty factor 
related to the unknown moment of inertia 2 ), and it does 
not depend on any details of the sparking gap model. It re- 
flects the fact that both the subpulse drifting rate (due to 
E x B plasma drift) and the polar cap heating rate (due to 
back-flow bombardment) are determined by the same phys- 
ical quantity, which is the potential drop across the inner 
acceleration region just above the polar cap. 

The microscopic properties of PSG model require a 
more sophisticated analysis, like the one presented in our 
Paper II. Here we can give simplified but more intu- 
itive estimate of the screening factor 77 = {a/2)(P / P4) = 
(1/2tt)/(P/P 3 ) and the number of sparks N = (P4/P3) = 
na, using arguments based on the complexity parameter 
a — r p /h presented in the paragraph below equation (1). 



3 OBSERVATIONAL VERIFICATION 

Table 1 presents the observational data and the predicted 
values of a number of quantities for four pulsars, which we 
believe to show clear evidence of thermal X-ray emission 
from the spark-heated polar caps as well as they have known 
values of the circulational subpulse drifting periodicity. The 
predicted values of P4 and/or L x are computed from equa- 
tion (3). Errors in L x is taken from the observational papers 
or derived from the distance uncertainty (taken from Cordes 
& Lazio (2002), except the case of B0656+14 for which it was 
obtained by Brisken et al. (2003) using the pulsar parallax), 
whichever is greater. The relationship expressed by equa- 
tion (3) is represented by the solid curve in Fig. 1, with two 

2 Considering general relativity, the moment of inertia of a neu- 
tron star can be written as I = 0.21MR 2 /(l - 2GM/(c 2 R)), 
where M and R is the neutron star mass and radius, respectively 
(Ravenhall & Pethick, 1994). Taking M = 1.4 solar masses and 
R ranging from 8x 10 5 to 1.7 x 10 6 cm, for the softest and stiffcst 
equations of state, respectively, one obtains the moment of inertia 
ranging from 7.82 X 10 44 to 2.25 X 10 45 g cm 2 , respectively. 
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Table 1. The observed and the model parameters for the four pulsars. 
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Note: t As P3 was not measured for these two pulsars we used the estimate of r] to calculate P3 /P. 



dashed curves describing the uncertainty in determining of 
the neutron star moment of inertia 2 . To save space in Table 
1 we give the basic pulsar parameters (P, P_ 15 , E x 1CT 32 , D) 
next to the pulsar name in the paragraphs describing each 
case below. 

PSR B0943+10 (P = 1.099 s, P_i 5 = 3.49, E = 
1.04 x 10 32 erg/s, D = 0.631 ^4 kps) is the best stud- 
ied drifting subpulse radio pulsar. As this case, along with 
PSR B1133+16, was discussed earlier in Papers I and II, we 
do not find it necessary to review it again (see Table 1 and 
Fig. 1). Error bars for P4 were given by Rankin & Suley- 
manova (2006), while errors for L x were given by ZSP05. 
See section 4 for discussion. 

PSR B1133+16 (P = 1.188 s, P_i 5 = 3.73, E = 
0.88 x 10 32 erg/s, D = 0.35±g;g| kps) is almost a twin of 
PSR B0943+10, in both radio and X-ray bands as it was 
demonstrated in Papers I and II, (see Table 1 and Fig. 1). 
Error bars for P4 were given by WES06, while errors for L x 
were given by KPG06. See section 4 for discussion. 

PSR B0656+14 (P = 0.385 s, P_ 15 = 55.0, E = 
381 x 10 32 erg/s, D = 0.288lg;gp kps) is one of the fa- 
mous Three Musketeers, in which the thermal X-ray emis- 
sion from the hot polar cap was clearly detected (De Luca et 
al. 2005). This pulsar is very bright, so the photon statistics 
are good enough to allow identification of the BB component 
in the spectrum. As indicated in Table 1, the X-ray lumi- 
nosity of this hot-spot BB component is L x ~ 5.7 x 10 31 
ergs/s. This value, when inserted into equation (2), returns 
the predicted value of P4 = 20. 6P. Amazingly, Weltevrede 
ct al. (2006b) reported recently the periodicity of (20 ± 1)P 
associated with the quasi-periodic amplitude modulation of 
erratic and strong emission from this pulsar. Thus, it is 
tempting to interpreted this period as the circulation time 
P4. Since there is no doubt about the thermal polar cap 
emission component, this case greatly strengthens our ar- 
guments given for PSRs B0943+10 and B1133+16, and the 
equation (3) receives a spectacular confirmation. It is inter- 
esting to note that the erratic radio emission detected by 
Weltevrede et al. (2006b) is similar to the so-called Q-mode 
in PSR B0943+10. The low frequency feature in the fluctu- 
ation spectra, identical to the one in the organized B-mode, 
was found by Rankin & Suleymanova (2006; see their Fig. 6). 
Asgekar & Deshpande (2001;AD01 hereafter) also detected 
this feature in the 35-MHz observations of PSR B0943+10 
(see their Figs.l and 2). This simply means that the E x B 
plasma drift is maintained in both regular (with drifting sub- 
pulses observed) and erratic (no drifting subpulses) pulsar 
emission modes. See some additional discussion in section 4. 
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Figure 1. The efficiency of thermal X-ray emission from hot po- 
lar cap L x versus circulation period P4 of drifting subpulses in the 
radio band. The solid curve represents the prediction of the PSG 
model (equation 3) with I45 = 1, while the dotted curves corre- 
spond to uncertainties in determining of the moment of inertia 1 . 
Error bars on P4 and L x were given by the authors mentioned in 
the text. 



PSR B0628-28 (P = 1.244 s, P_i 5 = 7.12, E = 
1.46 x 10 32 erg/s, D = 1.444l°; 2 7 ? kps) is an exceptional 
pulsar according to TO05. Its X-ray luminosity exceeds the 
maximum efficiency line derived by Possenti et al. (2002) by 
a large factor. However, one should note that PSR B0943+10 
with its luminosity derived from the PL fit, also exceeds 
this maximum efficiency (ZSP05, TO05). The BB efficiency 
L x /E ~ 1.9xl0~ 2 gives the predicted value of P 4 ~ (6±1)P 
from equation (2). It is very interesting that WES06 re- 
port the periodicity of (7 ± 1)P (see Table 1 and Fig. 1). 
According to the model expressed by equation (3) this rel- 
atively low modulation periodicity (i.e. high modulational 
frequency) can be interpreted as the circulation time P4. If 
this is true then PSR B0628-28 is not an exceptional pulsar 
at all. It lies on the theoretical curve in Fig. 1 at exactly the 
right place. This also means that the observed drift is highly 
aliased in this pulsar, with P3/P being considerably lower 
than 2. As concluded by WES06, this might be the case for 
most pulsars. Therefore, all or most features in modulation 
spectrum frequencies below about 0.2 cycle/P may in fact 
represent directly the E x B plasma circulation around the 
pole rather than the apparent subpulse drift periodicity. 
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4 CONCLUSIONS AND DISCUSSION 

Within the partially screened gap (PSG) model of the in- 
ner acceleration region in pulsars developed by GMG03, we 
derived a simple relationship between the X-ray luminosity 
L x from the polar cap heated by sparks and the circulation 
time Pi of the spark-associated drift detected in radio band, 
not necessarily in the form of regularly drifting subpulses. 
This relationship expresses the fact that both E x B drifting 
rate and polar cap heating rate are determined by the same 
value of the available potential drop. In PSRs B0943+10, 
B1133+16, B0628-20 and B0654+14, which are the only 
pulsars for which both L x and Pa are known at the moment, 
the predicted relationship between observational quantities 
holds very well (Fig. 1 and Table 1). This suggests that the 
PSG model may indeed be a reasonable description of the in- 
ner accelerator near the polar cap region. With the abundant 
radio drifting data (WES06) and the growing number of old 
pulsars detected in X-rays by XMM-Newton and Chandra, 
the clean prediction from this model (equation 3) will be 
unambiguously further tested with more pulsars in the fu- 
ture. PSR B0826-34 with P 4 about 14 or 7.5 P (Gupta, Gil, 
Kijak et al. 2004) and PSR B0834+06 with P 4 about 15 P, 
will be examined in the near future. 

For the carousel circulation time P4 to be measurable at 
all, it requires a strong unevenness in the circulating system, 
maybe a distinguished group of adjacent sparks or even just 
a single spark (see also scenario discussed by Gil & Sendyk 
(2003). Moreover, it requires this feature to persist much 
longer than the circulation time. Such favorable conditions 
do not occur frequently in pulsars and therefore direct mea- 
surements of P4 are very rare. In principle, in a clean case, 
using the fluctuation spectra analysis, one should be able to 
detect the primary feature P3, reflecting the phase modu- 
lation of regularly drifting subpulses, flanked by two sym- 
metrical features corresponding to slower amplitude modu- 
lation associated with carousel circulation. PSR B0943+10 
was the first pulsar to show such a model behavior (DR99) 
and PSR B0834+06 was the second one, as demonstrated 
by Asgekar & Deshpande (2005). The latter authors have 
also found a direct long period circulational features in both 
pulsars. For the B0943+10 they found it in their 35-MHz 
observations (AD01). In the case of B0834+06 Asgekar & 
Deshpande (2005) found an occasional sequence of 64 pulses 
with much weaker frequency modulation (present in the rest 
of their data) but with strong long period feature associated 
with the amplitude modulation due to the circulation of one 
or few sparks (see their Fig. 3)0 Most interestingly, how- 
ever, Rankin & Suleymanova (2006) were able to detect a 
long period circulational feature P4 in the so called Q-mode 
erratic emission mode in B0943+10. This apparently first 
detection of the Q-mode circulation time is very important. 
Indeed, this fact and other cases discussed in this paragraph, 
strongly suggest that no matter the degree of the organiza- 
tion of spark plasma filaments at the polar cap, the E x B 



3 A close inspection of this sequence of 64 pulses shows also a 
presence of even-odd modulation corresponding to the value of 
P3/2 close to 2. However, the slope of the secondary drift-bands 
changes sign, meaning that P3/P oscillates around the value of 
2 every P4 periods. Thus, at least in this sequence the subpulse 
drift in PSR B0834+06 seems aliased. 



drift motion is always performed at the same rate in a given 
pulsar. The problems is how to reveal this motion. 

Different methods of analysis of pulsar intensity fluc- 
tuations are sensitive to different effects. The method used 
WES06 has an obvious advantage of finding periodicities 
even in a very weak pulsars, so it resulted in a large increase 
of pulsars with drifting subpulses and/or periodic intensity 
modulation. Generally, WES06 can find only one period and 
they denote all the periods they find by P3, suggesting that 
these are primary drift periodicities. It does not have to 
be this way at all. In fact, we suggest that at least in three 
cases their reported values correspond to carousel circulation 
times P4,. We base our argument mainly on the fact that they 
satisfy nicely our empirical relationship (Eq. 3 and Fig. 1), 
without any obvious selection effect involved. Moreover, in 
B1133+16 the value of P 4 = (33±3)P is close to (37.4±?;*)P 
detected in the twin pulsar B0943+10 (see Paper I for more 
detailed discussion). In B0656+14 the periodicity of about 
20P results from intensity modulation of erratic spiky emis- 
sion, similar to the case of Q-mode in B0943+10. 

Using a concept of the complexity parameter a (GS00) 
corresponding to the ratio of the polar cap size to the spark 
characteristic dimension, we estimated a number of sparks 
,N, operating in the inner accelerating regions, as well as 
values of the screening parameter rj for the pulsars discussed 
in this paper. In the two twin pulsars both N and 77 are 
almost the same. Is seems trivial since in the approximation 
we used a depends only on the P and P values, which are 
close to each other for these two pulsars. However, N can 
also be found from the ratio of observed values of P4 and 
P3, and both estimates are consistent with each other. PSR 
B0628-28 seems quite similar to the twin pulsars, while in 
PSR B0656+14 the number of sparks is 4 times greater, 
and the screening parameter is quite high (corresponding to 
about 75 % of the vacuum potential drop). This is a result 
of relatively low P (large polar cap) and unusually high P. 
Thus, either the actual number of sparks is really that big 
in this pulsar, or the approximation of GS00 is not good for 
such a non-typical pulsar. It is not difficult to lower the value 
of the complexity parameter and a corresponding number of 
sparks (N = na) by a factor of 2 — 3, by considering larger 
radii of curvature of the actual surface magnetic field lines, 
or even the inverse Compton scattering instead of curvature 
radiation as seed photons for the sparking discharges (see 
Zhang, Harding & Muslimov 2000 and Gil & Melikidze 2002 
for some details). 
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